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Abstract 

5h , 

| Physical (and weak) regularity conditions are used to determine and classify all the possible types of 

spherically symmetric dust spacetimes in general relativity. This work unifies and completes various earlier 
results. The junction conditions are described for general non-comoving (and non-null) surfaces, and the 
limits of kinematical quantities are given on all comoving surfaces where there is Darmois matching. We 
show that an inhomogeneous generalisation of the Kantowski-Sachs metric may be joined to the Lemaitre- 
Tolman-Bondi metric. All the possible spacetimes are explicitly divided into four groups according to 
topology, including a group in which the spatial sections have the topology of a 3-torus. The recollapse 
' conjecture (for these spacetimes) follows naturally in this approach. 

o : 

O I. Introduction 

oo : 

0^ ■ The convenient Lemaitre-Tolman-Bondi (LTB) exact solutions have been exploited as the main inhomogeneous 

models in relativity and cosmology for many years. The remarkably rich structure of these solutions has many 
subtleties, in particular concerning the regularity of the metric ^, || ||. The purpose of this paper is to clarify, 
$1^ | unify and complete existing results on regularity. This topic has important implications, for example in the 

modelling of voids and other local inhomogeneities (see, e.g. Q). The new results include the matching between 
, an exact solution in the generalised Kantowski-Sachs family and the LTB solutions, and junction conditions for 

arbitrary non-null surfaces. The possible types of centre are extended. (Note that SS dust solutions need not 
possess a centre JtJ.) All possible composite spherically symmetric (SS) dust models are found and classified into 
four topological groups. One such topology, in which spatial sections have the topology of a 3-torus, appears 
not to have been discussed previously. 

We use comoving spatial coordinates, since these are best adapted to matching problems. For an analysis 
of astronomical observations, coordinates based on the past light cones of the observer are a better choice |l4| . 
The paper is laid out as follows. In section II, all of the solutions of Einstein's field equations (with a smooth 



SS dust source) are expressed in the forms (10) and (|14[). The geometrical requirements at any junction in 
spacetimes composed from these solutions are analysed in section III. In section IV, the reasonable physical 
regularity conditions are made explicit. This leads to regularity conditions within the domain of any one solution 
(section V), at any centre of symmetry (section VI), and at any junction between solution domains (section VII). 
Also in section VII, all the possible composite models are listed. In section VIII, the results are summarised 
and used to give a simple proof of the recollapse conjecture for these models. 

II. Spherical Dust Solutions 

In standard comoving coordinates x a = {t, r, 8, <p}, the dust 4- velocity is u a = 8f and the metric may be written 
as | 

ds 2 = -dt 2 + X(r, t) 2 dr 2 + R(r, t) 2 dn 2 , (1) 
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where dfl 2 = dO 2 + sin 2 9 dcj) 2 and we choose R > 0. The energy-momentum tensor is T b = /J,u a u b = —fiS^, 
where \i is the proper matter energy density. The Einstein tensor is j^] 
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where an overdot denotes u a d a — d/dt and a prime denotes d/dr. The units satisfy G = c = 1. 

In order for solutions of Einstein's equations G b a = 8irT b to exist, the Einstein tensor G h a must be defined 
through (^|)-(||). In such regions of spacetime, it is therefore required that 



X R; X and R are C 2 in f; X and R are C 1 in r; R is C 2 in r. 



(G) 



However, these regions may be joined together to form a composite spacetime, in which G\ is not defined by 
(H||) on the boundaries. 

Bondi integrated the system as follows. The G\ = field equation integrates to (assuming R' ^ 0) 



X = 



R! 



where E is an arbitrary function. The remaining independent equations reduce to 
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with M another arbitrary function. The corresponding proper density is given by 

M' 



M = 
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There are five solutions of (\8h 
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[0 < 77 < 2?r], 



with T(r) a third arbitrary function, t = et — T, e = ±1, and we denote the five solutions by (sl),..,(s5). 
No physical solutions exist for (E < 0, M < 0) or for (E = 0, M < 0). Note that (si) and (s2) are (locally) 
Minkowski spacetime. Motivated by equation (||), Bondi (|] describes M as a relativistic generalisation of 
Newtonian mass, and \E as a total energy. The surfaces t — are spacelike singularities. 

The above integration necessarily requires R' 7^ 0. If instead R' = 0, then a different solution results, as 
follows. The G r r — equation integrates to 
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where M > is an arbitrary constant and R =/= has been assumed, since the converse immediately leads to an 
inconsistency. The second integration (assuming < R < 2M, otherwise R 2 < 0) reveals that 

(12) 



(13) 



R = M (1 - cos 77) , 77 - sin 77 = Ar 1 Ut - Tj , 

where < 77 < 2tt and f is an arbitrary constant. By (|l2|), equation (^) implies the linear equation 

d 2 X X 



dvj 2 COS 7] — 1 



= 0, 



which transforms to a first-order Ricatti equation under X — ► X d n X. Hence the general solution of ( |l3| ) 
may be found provided one solution is known One particular solution is X = sin 77 (1 — cost;)" 1 , and the 
general solution is: 



(s6) 
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where A and B are arbitrary functions, and we denote this solution by (s6). Finally [by (|^)] the density reduces 
to 
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This solution is an SS variant of an inhomogeneous generalisation of the k = +1 Kantowski-Sachs metric, as 
was discovered previously [|| . The form of the solution presented in jl2| (in which \i ^ was assumed) is 
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The coordinate transformation between ( |l6| ) and (s6) is given by t = R, r' = A/C, f = (with the identifications 
C = 2A/B and a = 2M). 

In summary, the six possible SS dust solutions are (sl)-(s6). 



III. Junction Conditions 

The differentiability conditions (|^) on the metric need not be satisfied at the interfaces between domains of 
separate solutions in a composite model. Instead (weaker) matching conditions (for the geometry) must be 
satisfied |lC| . At these junctions, it is assumed that the coordinate r remains regular, which is not the case at 
a shell-crossing (by definition). Such singularities are excluded in this paper. The matching of two general SS 
spacetimes has been considered recently in Q , in which necessary conditions for the matching were presented, 
which are valid for any equation of state. Here necessary and sufficient conditions are found in the special case 
of dust. 

Consider firstly the case of a comoving boundary r = constant (which was analysed previously ||). The 
unit outward normal is n a — \X\S^, the metric intrinsic to the surface is 

Kb = 9ab - n a n b = diag (-1, 0, R 2 , R 2 sin 2 6) , (17) 

and the extrinsic curvature is 

K ab = h c a h d b X7 d n c = diag ^0, 0,^,^ sin 2 e) . (18) 

The Darmois matching conditions j^, ^) state that h ab must match across the surface and any discontinuity 
in K a b gives rise to a surface-density layer as represented by the Lanczos equation fl3l 

8^ ( 3 T ab ) = g cd (AK cd ) h ab - AK ab , (19) 
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where A denotes the limit of a quantity on the 'r+' side of the interface, minus the limit on the 'r— ' side. The 
surface energy-momentum tensor 3 T ab can be expressed in perfect fluid form 3 T ab = ( 3 /i 
and 3 p the surface energy density and surface (isotropic) pressure. Now 



V) u a u b + 3 ph ab with 



V = 3 T ab u a u b = —A ( 7 ^ T 
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Following Bonnor |J, this equation of state is regarded as unphysical, i.e. we require the extrinsic curvature on 
comoving surfaces to be continuous in r. By p7| ) and (|l8|), the junction conditions reduce to 

(20) 
(21) 



R continuous in r, 
R' 



\X\ 



continuous in r. 



For solution (s6), R! = 0, hence (s6) may only be matched (across a comoving surface) to one other SS dust 
solution, i.e. (s5) [by (g), {R! ->• 0, A ^ 0} requires J5 -> -1]. From (jl0|), (|l§) and @ the matching also 
requires M —> M and T — * T in the (s5) region. This motivates a characterisation of (s6): solution (s6) may be 
characterised within the LTB family by the conditions {M' = T' = 0, M > 0, = — 1}. 

Consider now the junction conditions on the spacelike surfaces t = constant. In this case, the unit normal 
is u a and the intrinsic and extrinsic curvatures are 

Kb = gab + u a u b = diag (0, X 2 , R 2 ,R 2 sin 2 0) , K ab = h c a h d b V d u c = diag (fj, XX, RR, RR 



(22) 

Hence continuity of h ab and K ab in t merely implies that R, X, R and X are continuous in t. (Spacelike surface 
layers, which imply an instantaneous transition, are treated as unphysical a priori.) In fact, by inspection of 
equations (|7|), (g) and (pi]), the metric tensors for (sl)-(s6) are infinitely differentiable (C°°) in t. 

An analysis of the Darmois matching problem for solutions (sl)-(s6) across a general non-comoving (and 
non-null) SS surface provides insight into the nature of general dust models. From (Q), the unit normal n a and 
unit SS tangent t a to such a surface satisfy 



-t a t a = A, 



i a = (-p, \x\Vx + p 2 , o, o) , 



ta 



where the surface is timelike (spacelike) for A — 1 (A = - 
the surface: 

, . ds 



VX + P 2 , P\X\, 0, Oj , 
1) and P is a function determined by the equation of 
P 



dt \X\VX + P 2 
A lengthy calculation leads to the extrinsic curvature of the surface: 
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Since the coordinates continue through the surface, all four coordinates are induced on the surface, and one of 
r, t is redundant there. Now, the Darmois conditions with no surface layers take the form 



Ah ab = => AR = A\X\ = 0, 
AK ab = => AR' = A (x/x) = 0, 



(26) 
(27) 
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where it has been assumed that P ^ and P ^ 1, as these cases are already considered above. By equations 
(H), (H) and d|), AP = 0. We now show that 

If there are no surface layers, then all boundaries between domains of different solutions (sl)-(s6) in 
a composite dust model must be comoving, i.e. {r — constant}. 



(s5), \X\ = \R'\/VTTE. Then by 



The proof is as follows. In matching together two different solutions (si) 
AE = 0. From © 

k= Rgt a _ PR' 

y/X + P 2 \X\V\ + P 2 ' 

so that impiy AR = °- Thcn b Y ®, AM = °- Finally AR = forces AT = 0. Hence all three LTB 

functions E, M and T 'carry through' the surface. Since E, M and T depend only on r, the surface must be 
of the form {r = const}. By (p7|) , matching between one of (sl)-(s5) and solution (s6) requires R' — > on the 
(sl)-(s5) side. This is only possible in (s5), with E — > — 1. However [by (0)], E(r) = —1 only at isolated values 
of r, i.e only on a comoving surface. The result follows. 

It is of some interest however to establish the nature of non-comoving singular surfaces (surface layers) in 
these solutions, in which case Ah a b = ^ AK a b- Restricting to A = 1 (timelike surface layers) one may 



construct 3 T a b once again, by fll9|) 
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F, = 
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A comparison between ( p8| ) and the perfect fluid energy tensor 3 T a b = ( 3 M + a p) t a tb + 3 ph a b reveals that the 
surface layer is always of perfect fluid form. The surface energy density and pressure are 



AR' 
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4irR\X\y/l + P 2 ' 
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2AR' 



The nature of the matching problem changes greatly in moving away from the comoving case. In the non- 
comoving case, conditions at the wall must be satisfied through some range of r. These conditions are therefore 
'dragged' into the adjoining spacetimes, since the arbitrariness in these solutions resides purely in functions of r. 
This approach to surface layers and its application to models of voids will be the subject of a subsequent paper || . 

In the rest of the current paper, singular surfaces are ruled out. All boundaries are necessarily comoving, 
and the metrics are fully determined then by choices of the arbitrary functions E(r), M(r) and T(r) [and A(r) 
and B(r) in regions where E(r) = — 1]. Throughout these SS dust models, junction condition (|20|) reduces to 



M(r)(> 0), E(r)(> -1) and T(r) are continuous. 



(29) 



IV. Regularity Requirements 

In this section, the physical requirements (and one coordinate constraint) to be imposed on the metric are made 
explicit and justified. From section III, the metrics of (sl)-(s6) are C°° in t, hence attention is focused on radial 
differentiability. 

For a well-behaved radial coordinate, g rr must be piecewise continuous in r (that is, continuous except at 
isolated values of r, where both left and right limits must be finite). It is also required that lim± g rr 
everywhere, where a + (— ) denotes a right (left) limit. This extra condition purely defines a 'good' spatial 
coordinate, so that dr is everywhere proportional to the differential increase in radial proper distance. This 
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subtlety ensures that the continuity properties of physical quantities may be expressed unambiguously through 
their differentiabilities in 7'. 

The dust is characterised (in the SS case) by the density fi, expansion rate O = K a a and shear a a b — 
Kab — 3@^o6- For physically reasonable matter, fx, O and a a b must each be piecewise continuous in r. By 
spherical symmetry, a a i, — * wherever R — > 0. For the metric (fy) 



b 1 



X R 

X ~ R 



xdiag (0,2, -1,-1), ©-2| + ^, (30) 



where X/X and R/R are 'radial' and 'azimuthal' expansion rates. 

To ensure that the spacetime itself is regular at each point jl5j it is required that 
R(i)(j){k){i) — -^afccde^e^e^e^ is piecewise continuous in r, where e%<. is an orthonormal tetrad basis. Here 
a natural choice is made: e? % = u a , ef^ — \X\^ 1 S'^, e%-, — R~ 1 5@, e%-, — R~ 1 cosec95^. For solutions (sl)-(s5), 
the nontrivial components are found to be 

2M M 

«(0)(1)(0)(1) = 47T/J - — , i£( )(2)(0)(2) = -^3 , 

M „ 2M 



#(1)(2)(1)(2) = 47T/J- -^3, i?( 2 )(3)(2)(3) - ^3 



and for (s6) 



XR M _M 

«(0)(1)(0)(1) - YR ~ #3' K (0)(2)(0)(2) - -^3, 

XR _2M 

«(1)(2)(1)(2) - rt(2)(3)(2)(3) - -£§-• 

Hence, for regular spacetimes: at points where i? — > 0, lim(M/i? 3 ) must be finite [using (|29|), and since /z, X/X 
and R/R are already required to be piecewise continuous in r\. Solution (s6) does not admit central points. To 
summarise, throughout the models it is required that: 

Rl. The junction condition ( pl| ) is satisified, 

R2. g rr is piecewise continuous in r and lim± g rr =/= 0, 

R3. /i > is piecewise continuous in r, 

R4. R/R is piecewise continuous in r, 

R5. X/X is piecewise continuous in r, and X/X — > _R/i? wherever R — > 0, 

R6. M /i? 3 is finite wherever R —> 0, 

(except, trivially, at the spacelike singularities r — > 0). 

In sections V-VII the above conditions [with ([29]) satisfied a priori] are enforced for a general SS dust metric, 
to guarantee regularity. 

V. Regular Solutions 

Here the conditions of section IV are verified in turn, for points in the domain of a solution. This domain does 
not include the origin (treated in section VI) or interfaces between solutions (treated in section VII) . 



Rl. Junction condition ( pl| ) is automatically satisfied for (s6). For any of (sl)-(s5) it implies that [using (Q)] 
R' may change sign only at values of r satisfying E{r) = —1, (31) 



G 



as was noted in ||. Hence R' > throughout the domain of each solution (sl)-(s4), or R' < throughout. In 
(s5), R' may change sign where E = — 1. Now for (s4), i?' may be written as 
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M' 
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77 sinh 77 
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(32) 



and for a positive density [by equation (^)] , R' and M' must have the same sign. At large rj (large t) , the third 
term in ( |32| ) dominates, so that R' and E' must have the same sign. At small 77, the second term dominates, so 
T" must have the opposite sign to R'. These conditions are also sufficient (to ensure R' has a constant sign), 
since in (|3|) each function of 77 in square brackets is strictly positive for all allowed 77. Hence [for (s4)] equation 
( |3l| ) implies that 

± R' > {±M' > 0, ±E' > 0, ±T' < 0} . (33) 



For any of (sl)-(s3), similar reasoning also leads to (|33|). Now for regions of (s5) which satisfy E ^ —1 (in which 
R' cannot change sign), it is useful to write R' in the form 
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3E 
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+ ^[1- cos 77] 



(34) 



where each function in square brackets is always positive. At small 77 the second term dominates. As 7/ — > 2ir 
the first term dominates, and M' must have the same sign as R' for /1 > 0. Hence in these regions of (s5) 



±R' >0^> { ±M' > 0, ± E'M 



2 T'\E\ b / 2 
-M'E+ 11 

3 3tt 



> 0, ±T' < 



(35) 



Equations ( j33| ) and ( pq ) are the Hellaby and Lake no-shell-crossing conditions ||, derived here from the junction 
conditions. Violation of ( (33| ) or ( |35| ) would necessitate either a pathological choice of radial coordinate or true 
caustic formation. 



R2. From (g), piecewise continuity of the radial coordinate within the domain of one of (sl)-(s4) [and in regions 
of (s5) where E(r) ^ — 1] requires piecewise continuity of R' in r, for which it is necessary and sufficient that 
M'(r), T'(r) and E'(r) are piecewise continuous. Now lim± g rr ^ implies lim± R' ^ in these regions, which 
[from ( j33|) and (^B])] forces one of lim+ M', lim+ T' and lim+ i?' to be nonzero at each point (likewise for the 
left limits). At points in (s5) where E(r) = — 1, lim+ g rr is finite if and only if 

M' T' E' 

lim_i_ — , lim_i_ — , and limj_ — , are finite, and at least one is nonzero, (36) 

WTTe' VT+e VT+e 

as follows from (g), ( p4] ) and (35) (and analogously for lim_(7 rr ). One necessary consequence of ( p^ ) is that 
lim± R' = lim± M' — lim± T' = lim± £" = at these points. Finally, for (s6) the radial coordinate is well- 
behaved if A(r) and B(r) are piecewise continuous and if at least one of lim + A and lim_|_ B is nonzero [by JT4|), 
and likewise for lim_l. 



R3. The density is piecewise continuous in (sl)-(s5) as a consequence of the junction conditions, which is seen 
as follows. From ([)[), fx is at least continuous (in r) except at isolated points, since R' and M' have this property 
(see above) and R' = only at isolated points. Trivially fi = for (si) and (s2). For (s3)-(s5), if M' = 
in a finite region, then /i = there. Otherwise lim± /i are finite if and only if lim± (R' /M') ^ 0, but this is 
automatically satisfied in (s3)-(s5) [e.g. in (s4) 



Rl 1 
lim-t — > — 
± M' ~ E 



rj sinh 77 



cosh 77 — 1 



>0, 



by (|32|),(p3|)1. Hence lim± fi are finite, and /i is piecewise continuous in r. In fact, lim± /i = (with M' ^ 0) only 
if either lim± (E' /M') = +00 or lim± (T'/M') = —00. Note that none of these remarks require any modification 
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at points in (s5) with E = — 1. Now the density in (s6) [given by (|l5|)l vanishes if B = and 4^0 [and (s6) 
degenerates to part of the exterior Schwarzschild solution]. Otherwise \x is finite and positive at all times in (s6) 
if and only if 

lim±->- (37) 

A 7T 

which ensures no zeroes in the denominator of ([l5]). Piecewise continuity of fi in (s6) follows from the piecewise 
continuity of A(r) and of B(r). 

R4. Full continuity of R/R'mr is guaranteed throughout the SS dust spacetimes by (||) and (p9|). 



R5. For (si), X/X = 0. For (s2)-(s5) [by (@)-(|)] 



X _ R 
X~R 
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Hence piecewise continuity of the shear in (sl)-(s5) requires the further condition that \im±(E' / R 1 ) is finite 
wherever R' — > 0. This is trivially satisfied in (s3). In (s2), (s4) and (s5) it is also automatically satisfied since 
Km±(R! /\E'\) > [e.g. in (s4) 



R' M 
lira + — > — ^ 

± E' ~ E 2 



sinh?/ (sinh?7 — 3?/) 
2 (cosh 77 — 1) 
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by (B2|),(B3)]. Finally, solution (s6) has 
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X M(l-cosr)) 



\B (77 - sin 77) - A 



B{\ — cos 77 — \r\ sin 77) + A sin 77 



which is automatically piecewise continuous in r by ( |37j ) and by the piecewise continuity of A(r) and B(r). 



VI. Regular Centres 



In this section the possible types of origin (for which R — > 0) are determined by imposing the conditions of 
section IV. Only comoving origins are possible, and they may join only to solutions (sl)-(s5). All the results are 
given in table |lj, in which (i) derives from the condition R — > 0, (ii) derives from lim-t g rr 0, and (iii) forces 
the shear to vanish. In each case MjR? — > |7T/i, which is the Newtonian limit. Examples are provided; in each 
the origin lies at r = 0. The central behaviours of (s4) and (s5) generalise previous results. 
(INSERT TABLE [l] HERE) 

One central limit suffices to illustrate the arguments used to obtain table |l|. At an origin of (s4), suppose 
77 — * 00 . Then by (ftCj) , on approaching the origin 



2E 3 ' 2 t 



M 

and R -> forces E -> 0, M -> 0. Now 

R e X ei^E'E- 1 ' 2 + T'E- 1 Mt- 2 +M'E- 1 t- 1 ) 



(38) 



R t ' X [M'E- 1 log (E^M- 1 ) - T'E 1 / 2 + \E'E- x I 2 t\ ' 
so that vanishing shear requires either 
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or 



, :e i / 2 e'\ ,. (MV . 

lim = lim = 0, lim 



log 



E 3/2 \ jv E 3/2 



M 



M> 



= 1. 



However the latter case ( |40| ) is ruled out since the logarithmic term must diverge, by 
only possibility, and this reduces to 



liit) ( 1 : lim 



M' /E 3 ^ 2 
E 1 / 2 E> ° g \ M 



0. 



by consideration of (|33|) and (38). From (32) 



R -> log - T'£ 1/2 + — . 



so that [with (f4l|)], lim<? rT . 7^ requires lim 1 / 2 ) finite and nonzero. Finally 



47T/L4 



r 2 log 



£ 3 /2X E 3/2 Tt 2 E 1/2 E ' t 3 



M J M' 2M' 
which vanishes (again by the divergence of the logarithmic term). 



Hence 



(40) 
|) is the 

(41) 



(42) 



VII. Regular Interfaces 

In this section the conditions of section IV are considered on the comoving interfaces between domains of 
solutions (sl)-(s6) in a composite model. The solution domains are assumed to be regular (in the sense of 
section V) and this generally ensures regular interfaces. 

By (fU|), the sign of R! cannot change across these interfaces, since E ^ —1 on them [except on interfaces 
between (s5) and (s6), but R! = in (s6)]. Solution (si) may not be matched to any other solution, since R = 
in (si) (whereas the other solutions are cosmological) . Solution (s2) does not match to (s3), because M — > 
forces R — > in (s3). Also, (s2) does not match to (s5), since (E — > 0, M — > 0) forces R — > in (s5). From 
section III, (s6) only matches to (s5). There remain just five physical types of junction, given below. At each 
of the five, equation (^) ensures continuity of R. 



a. Matching (s2) to (s4) 

The (s2) side of this interface is unconstrained by the matching. On approaching the interface from (s4), M — * 0, 
E 7^ and M' > throughout some finite region [by the piecewise continuity of M' , and since M > in (s4)]. 
Hence by (pi]) and (p3|), R must increase in the direction (s2)— >(s4). Now rj obeys (p8|), so that in (s4) 



1 



*9r 



VTTe 



' M' , 




^log 


( E v) 







— T E 



I W U2 



E't 
2E 1 / 2 



Therefore r is a good coordinate if lim.( S 4) T', lim( S 4) E' and lim( S 4) M' log M are finite, and at least one of them 
is nonzero. 

On the (s4) side, the density reduces to p2]), and vanishes by the divergence of the logarithmic term. [Note 
that fi = in (s2).] 

On the (s4) side, since (M —> 0, M'logM finite) forces M' — > 0, we have 



X 
X 



whereas on the (s2) side 



_ 2ET' _ 2M' , M 



X 



if E' -> 0, 
otherwise, 



X T-2ET'/E n 
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and on both sides R/R — > e/r. Hence the shear is necessarily finite on both sides of the interface, as required. 



b. Matching (s3) to (s4) 

The (s3) side is unconstrained by the matching. On approaching the interface from (s4), E — > while M 7^ 0, 
so that 77 w (Gr/M) 1 / 3 ^ 1 / 2 — > 0, and £" > throughout some finite region [since E > in (s4)]. Hence by @) 
and (|33]), i? must increase in the direction (s3)— >(s4). On the (s4) side 



M' 



6AP 



— T I - 



V 3r 



£1 

40 



M 



-,1/3 



so that r is a good coordinate provided 

lim( S 4)M', liiri^T' and \vaxu4\E are finite, and at least one is nonzero. 
On both sides of the interface, the density reduces to 

f 



and on both sides, 



47T/1 



X 
X 



3 T 2 

2 ' 



ZMT't _r_ (drfH^^E' 



mow 



if M' -> 0, 
otherwise, 



1 /3 

1 MV , / 243r 5 N 7 p/ 
M I 250M 2 J ^ 



3M't 2 



2M 



3T'r 
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R 
R 



2e 
3r' 



Therefore both fj, and X/X are well-behaved. 



(43) 



(44) 



(45) 



(46) 



c. Matching (s3) to (s5) 

This interface is similar to b. On approaching the interface in (s5), r\ (Br/M) 1 / 3 ^! 1 / 2 — > 0. Then all the 
results in b for the kinematics and radial coordinate (43)-(|46|) follow [with '(s4)' replaced by '(s5)']. 

However, R must increase in the opposite sense to that in b - here R must increase in the direction (s5)— *(s3), 
as we now show. Taking the limit of equation ( |35| ) on the (s5) side gives sign(i?')£" > 0. Since E < in (s5) 
(and E' is piecewise continuous), the result follows. Note that the coordinate condition (Eij) is crucial to this 
proof. Note also that since R/R — > > at any interface between (s3) and (s5), the r-continuity of R/R 

forces the existence of a finite region in (s5) adjoining the interface for which the azimuthal expansion rate is 
positive R/R > [even though all points in (s5) eventually satisfy R < 0]. 



d. Matching (s4) to (s5) 

Both sides are constrained by E — > with M > 0, and the resulting junction is given precisely by combining 
the results of b and c. In this case, R must increase in the direction (s5)— >(s4), by a similar argument to that 
given in b. 

e. Matching (s5) to (s6) 

At this interface the (s6) side is unconstrained by the matching. Equation (|3l]) places no restriction, and R may 
increase in either direction on approaching the interface from (s5). Now lim( s5 ) E = — 1. Therefore, as discussed 
in section V, r is a good coordinate provided ( |36| ) is satisfied. 

The (s5)-limits of /j, and X/X are just those of an 'ordinary' point, i.e. one in the domain of (s5). In this 
sense, the matching conditions at this type of interface are considerably less restrictive than those at the other 
four. 
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Combining these results with the rest of the paper, all the regular SS dust models may be classified into four 
topologies: 



i. Open models with one origin 

By noting the sense in which R must increase at the interfaces a-e above, the only possible composite models 
are: 

G(sl) + , 0(s2) + (s4) + , 

G(s2) + , 0(s3) + (s4) + , 

C(s3) + , C(s5) + 5(s5) + (s3) + , 

C(s4) + , C(s5) + 5(s5) + (s4) + , 

C(s5) + 5, C(s5) + 5(s5) + (s3) + (s4) + , 

where O denotes an origin, and a superscript + (— ) implies that R increases (decreases) from left to right. 
Here S is any combination of (s5)~, (s5) + and (s6). Note that open models can be constructed from collapsing 
solutions [e.g. 0(s5) + (s6)]. Papapetrou jl6j discussed a particular example of C(s5) + (s3). 

In the above construction, we have noted from (Q) that on t =const, d\ = \dR\/y/\ + E, where \ is radial 
proper distance. Hence by (pl|), if E > a > —1 for all \ > P ( a ,f3 constants) then: 

dR 

X —> oo forces R — » oo if — > 0, (47) 

dx 

dR 

there is a finite value of x > P f° r which R = 0, if — < 0. (48) 

dx 

However, if E — ► — 1 as x oo, then neither of (f47j) , (fl8|) are necessary. An example of C(s5) + (s5) of type i is 



sin r 
7*1 



E={ s^Voj. 1 -^ ] fOT 0<r<r , M= ^[M,+e^] for < r < r , 
) — 1 + e 2r for r > r , \ M^ + e r for r > r , 

T = 0, Moo > 0, 7T < r < 2tt, (49) 



and 



E= [ [l-e" 2 ^] for 0<r<r o , M = \ ^ [M^ - e"^] for < r < r , 

I — 1 + e _2r for r > r , ) Afoo — e _r for r > r , 

T = 0, < Moo < 2/3, r > 0, (50) 

is an example of C(s5) + . In each of ( |49|) and (^o|), i? — >const> as x - * oo. There are no SS dust models with 
R — > as x — > 00 [by (0) and since, by (]l0|), R — > requires E 1 — > 0]. 



ii. Open models with no origin 

By (^8|), to avoid a zero in R, a model with no origin must either be composed entirely of (s6), or it must contain 
a section of (s5), in order to allow (at least one) minimum in R. Then the possible matchings are evident: 

(s3)-(s5)- 1 J (s5) + (s3) + 
(s4)-(s5)- r) (s5) + (s4) + 
(s4)-(s3)-(s5)- J I (s5) + (s3) + (s4)+ 

Examples and a detailed analysis of such models are provided in 0. In these models, due to the presence of 
collapsing solutions (s5),(s6), an origin does eventually form, but gravitational collapse will violate the regularity 
conditions in any case. 

iii. Closed models with two origins 

These models must contain a section of (s5), since there must be (at least one) turning point in R. The models 
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cannot contain a section of (s2), (s4) or (s2)(s4), since the section would either contain an origin and match to 
another solution, or would match to other solutions on both sides. Hence E would vanish on both sides, and 
since E > throughout the domains of (s2) and (s4), E' could not have the same sign throughout, contrary to 
)]. There can be no (si) section in the closed model, since it does not match to any other solution. 
There can be no (s3) region in the model either, since R must increase in the direction (s5)— »(s3). Hence if (s3) 
contains an origin, it cannot match to (s5). Conversely, if (s3) does not contain an origin, it cannot match to 
(s5) on both sides, leaving the model open. This leaves just (s5) and (s6) to construct these models, and the 
possibilities are: 

0(s5) + S(s5r0 



(H) [with (31 



iv. Closed models with no origin 

Consider an SS dust model which has R > in some range < r < d (and at some t). This final possibility of 
composite models is obtained by identifying (matching) the surfaces r = and r = d. Since AR = 0, the model 
must be everywhere (s6) or else it must contain a section of (s5) [otherwise sign(i?') is constant in < r < d, 
which forces i?(0) ^ R(d)]. No sections composed from the solutions (sl)-(s4) may be present, since they would 
be forced to match to (s5) on both sides. This would force R' to change sign in the section (since R must 
increase away from (s5) into these solutions) and this is not possible, by (|3l"|). Hence the models may only be 
constructed from (s5) and (s6), with the possibilities: 

1ST 



where I denotes the surfaces which are identified (at which the standard matching conditions must be satisfied, 
as we have described). The spatial sections of these models have the topology of a 3-torus. An example is 
provided by 



ar 



1 



for < r < id, 



E = < 



a(r - \d) 2 - 1 for id < r < |<2, 



b + cr 2 for < r < jd, 

M=\ b+ led 2 - c(r - irf) 2 for \d < r < §d, 



^ a{r~d) 2 



1 for jd < r < d, 



b + c(r - df 



for jd < r < d, 



1 



T = 0, a 2b+-cd 2 < -c 



(51) 



where a,..,d are positive constants. Note that a model of type iv cannot be constructed from the homogeneous 
(Friedmann-Lemaitre-Roberston- Walker) subclass of LTB (since the elliptic homogeneous solution has only one 
point with E = —1, at which R is maximum). 



There are no further possible topologies or composite models. Examples of models of types i-iii are given in 
previous literature (see especially Q). 

VIII. Conclusion 

In this paper, full regularity conditions have been derived and discussed for SS dust spacetimes. From section 
HI, the solutions for the metric in section II (which are all C°° in t) may be joined only on comoving surfaces. 
Hence the models are fully determined by choices of the functions E(r) > —1, M(r) > 0, T(r), A(r) and B{r). 
Existence of the Einstein tensor places a basic restriction on these functions [equation (Q)]: E, M and T must 
be C 2 , A and B must be C , except at a finite number of points. At these points, there may be discontinuities 
in any of E' , M', T', A or B. However there are additional, more subtle requirements of the functions. For 
example, at points where E —* —1, equation ( ]36| ) must hold. These differentiability conditions ensure the good 
behaviour of all the relevant physical quantities. 

In the current context, the recollapse conjecture [O, i.e. 'all closed SS dust models must recollapse ev- 
erywhere', follows directly from the results of section VII. We simply note that all the possible closed models 
contain only the solutions (s5) and (s6), which both recollapse in a finite time. This is a simple alternative 
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proof to that of Burnett g (which involved considering the lengths of timelike curves in these spacetimes). Our 
proof slightly strengthens that of Bonnor || , in that no mathematical assumptions are required (as were used 
in ||) other than those explicitly demanded by the regularity. 

Acknowledgements: We thank Bill Bonnor and Charles Hellaby for helpful discussions. 
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Soln. 


Behaviour of E, M and T 


Kinematics 


Example 


(si) 


(i) T —> 

(ii) limT" finite, nonzero 


6 = 
/i = 


T=-r 


(s2) 


(i) E - 

(ii) lim(E 1 / 2 E') finite, nonzero 

(iii) ET' jE' -> 


e -» 36T- 1 

/i = 


£ = r 2 
T = 


(s3) 


(i) M -» 

(ii) lim(Af" 2 / 3 M') finite, nonzero 

(iii) MT'/M' -> 


e -» 26T- 1 

47T/J — » |t~ 2 


M = r 3 
T = 


(s4) 


(i) E z l 2 /M -» 0, M -» 

(ii) lim(M" 2 / 3 Af) finite, nonzero 

(iii) lim(MT'/M') = lim(M 1 /3^7M / ) = 


e -» 26T- 1 

47T/J — > §t~ 2 


E = r 3 
M — r 3 
T = 


(i) E z ' 2 /M -» +oo, £ -> 

(ii) lim(i? _1 / 2 i?') finite, nonzero 

(iii) Um^T'/E') 

= lim [S- 1 / 2 ^/'^'- 1 log(S 3 / 2 M- 1 )] = 


e -» 36T- 1 

4717* — » 


£ = r 2 
M = r 4 
T = 


(i) E 3 / 2 /M — > a > 0, M — > 

(ii) lim(_E M') finite, nonzero 

(iii) E'^ME'/M' -» 2/3 and MT'/M' -» 


— > 3easinhr7/(cosh?7 — l) 2 
47T/J — > 3a 2 /(coshr? - l) 3 


£ = r 2 
M = r 3 
T = 


(s5) 


(i) |£f/ 2 /M — > 0, M — > 

(ii) lim(M~ 2 / 3 M') finite, nonzero 

(iii) lim(MT'/Af') = lim(M 1 / 3 £;7M') = 


6 -» 26T- 1 

47T/J — > §T~ 2 


£ = -r 3 
M = r 3 
T = 


(i) |£| 3 / 2 /M -» a > 0, M -> 

(ii) lim(i? _1 M') finite, nonzero 

(iii) E~ X ME' /M' -> 2/3 and MT'/M' -» 


6 — > 3ea sin ?]/(l — cosi]) 2 
Air/I — > 3a 2 /(l - cost;) 3 


£= -r 2 
M = r 3 
T = 



Table 1: Central Behaviour 
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